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Abstract:  A  computational  model,  based  on  diffusion,  differential 
rotation,  and  meridional  flow,  has  been  developed  to  simulate  the 
tranport  of  magnetic  flux  on  the  Sun.  Using  Kltt  Peak  magneto¬ 
grams  as  Input,  we  have  determined  a  best-fit  diffusion  constant  by 
comparing  the  computed  and  observed  fields  at  later  times.  Our  value 
of  730  +  250  km^/s  Is  consistent  with  Leighton's  (1964)  ^tlmate 
of  770-1540  km2/s  and  is  significantly  larger  than  Mofh6r's  (1977) 
estimate  of  200-400  km^/s.  This  suggests  that  diffusion  may  be  fast 
enough  to  account  for  the  observed  polar^maghetlc  field  reversal 
without  requiring  a  significant  assist  from  meridional  currents. 

-^Thls  paper  presents  the  initial  results  of  a  project  to  simulate 
the  transport  of  solar  magnetic  flux  using  diffusion,  differential 
rotation,  and  meridional  flow.  Hie  study  concerns  the  evolution  of 
large-scale  fields  on  a  time  scale  of  weeks  to  years  ,  and  Ignores 
the  rapid  changes  that  accompany  the  emergence  of  new  magnetic  regions 
and  the  day-to-day  changes  of  the  supergranular  network  itelf. 

Our  Initial  objective  was  to  determine  the  value  of  an  effective 
diffusion  constant  that  would  provide  the  best  fit  between  the 
computed  and  observed  fields.  To  our  knowledge,  no  fully  quantitative 
determination  has  been  attempted  using  modern  computational  techniques 
and  high-quality  observations.  A  resolution  of  the  discrepancy 
between  Leighton's  (1964)  model-dependent  estimate  of  (770-1540 
km^/sec)  and  Mosher's  (1977)  semi-quantltatlve  estimate  of  (200-400 


km^/sec)  would  help  to  determine  the  relative  Importance  of  diffusion 
and  meridional  flow  in  transporting  flux  to  the  Sun's  poles.  With  a 
value  as  low  as  Mosher’s,  diffusion  vrould  require  a  substantial 
assist  from  meridional  flow  to  produce  a  timely  reversal  of  the 
polar  magnetic  field,  whereas  with  a  value  as  high  as  Leighton's, 
diffusion  could  accomplish  the  reversal  alone. 

We  assume  that  the  photospheric  field's  radial  component 
satisfies  a  continuity  equation  containing  the  transport  effects  of 
Leighton  diffusion,  differential  rotation,  and  meridional  flow.  In 
spherical  coordinates  this  equation  is: 

3B  1  f  1  3  ,  .  ,  33»  1  323,  ^  1  f  2a  3B, 

3t  ■  Ti  Uin9  30^^^"®  36^  sin^0  T2  ®  3(t>' 


[jjijij  (B  sine  sinSe) 


where  the  time  constants  areT^^  ~  R^/K  (K  is  the  diffusion  constant), 
T2  -  1/cJq  (Wq  is  the  differential  rotation  rate),  and  T3  ■ 

R/Vq  (Vo  Is  the  amplitude  of  Che  meridional  flow  speed) .  Typical 
values  are  Tj  ■  20  years  (for  K  ■  800  km^/sec),  TV  “  21  days  (for 
a  Newton  and  Nunn  rate  of  2.77  deg/day),  and  T3  ■  1  year  (for  a 
meridional  flow  amplitude  of  20  m/sec  (Dtivall ,  1977)). 


Our  study  has  progressed  through  several  Increasingly  difficult 
and  realistic  phases.  First,  to  test  our  modeling  procedures,  we 
used  photographic  prints  of  ICltc  Peak  magnetograms  to  estimate  the 
location,  pole  separation  and  strength  of  new  bipolar  magnetic  regions 
during  1976-1981  as  the  present  sunspot  cycle  evolved.  Using  these 
doublet  sources ,  and  an  asstmed  initial  global  dipole  field  of  1  Gauss , 
we  calculated  the  evolution  of  the  large-scale  fields  and  displayed  the 
fields  in  Carrington  synoptic  charts  which  we  updated  at  27-day  inter¬ 
vals  Co  include  the  contributions  of  flux  from  newly  emerging  regions. 

Second,  to  estimate  the  diffusion  constant,  we  selected  four 
relatively  isolated  magnetic  regions  which  were  observed  on  several 
consecutive  solar  rotations.  For  each  of  these  regions,  our  objective 
was  Co  deposit  the  KiCt  Peak  digital  magnetic  measurements  into  the 
simulator  during  the  Initial  solar  rotation,  and  to  compute  the 
evolving  large-scale  magnetic  fields  at  elapsed  times  of  1,  2,  and  3 
rotation  periods  for  a  range  of  values  of  the  diffusion  time  In 

practice ,  we  also  considered  other  combinations  such  as  depositing 
the  flux  during  its  second  disk  appearance  and  comparing  the  computed 
and  observed  fields  on  the  third  rotation.  This  procedure  not  only 
provided  additional  tests  of  our  measurement  consistency  but  also 
helped  to  avoid  possible  magnetic  field  calibration  errors  associated 
with  Che  presence  of  sunspots  during  the  initial  rotation.  In  this 
second  phase  of  our  study,  we  deposited  digital  data  for  only  a 
limited  area  including  the  source  region.  The  surrounding  area  con¬ 
tained  the  typical,  but  unrealistic,  background  field  computed  from 
the  doublet  sources  of  phase  one.  To  minimize  the  unknown  errors 


associated  with  this  unrealistic  background  field  we  limited  our 
comparisons  to  the  immediate  area  of  the  expanding  bipolar  region 


For  each  region,  the  best  fit  was  obtained  numerically  by 
minimizing  the  sum  of  the  squared  errors  between  the  ccMsputed  and 
observed  fields  and  by  allowing  for  a  possible  bias  In  the  zero  level 
of  the  Kltt  Peak  measurements.  A  relatively  broad  minimum  was  always 
obtained  probably  due  to  the  unknown  Influence  of  the  weak  background 
fields.  However ,  In  half  of  these  cases  the  minimum  solution  also 
provided  close  agreement  between  the  observed  and  computed  values  of 
the  centroids  and  the  total  surviving  flux  of  the  magnetic  region. 

The  results  for  the  case  of  VgK)  are  shown  in  the  following  table: 

TABLE  1  -  DIFFUSION  CONSTANT  MEASUREMENTS 
Region  ABBCCCDD 

Deposition  Rotation  10200102 

Comparison  Rotation  21312213 

Diffusion  Constant  470*  680*  700*  1160  1140  1220  1690  1070* 

Summary:  K  >  1050  +  400  km^/s  (all  8  measurements) 

K  a  730  +  250  kffl^/s  (4  best  measurements) 

*  most  reliable  determinations ,  based  on  considerations  of  the  loca¬ 
tion  and  total  flux  of  the  evolving  region. 

The  eight  determinations  in  Table  1  yield  an  average  diffusion 
constant  of  1050  +  400  km^sec,  but  the  four  most  reliable  measure¬ 
ments  give  the  smaller  value  of  730  +  250  km^/sec.  Both  are  signi¬ 
ficantly  larger  than  Mosher's  (1977)  200-400  km^/sec,  and  are  consis¬ 
tent  with  Leighton's  (1964)  770-1540  km^/sec.  In  future  phases,  we 
plan  to  study  the  effects  of  including  the  observed  background  field 
as  well  as  a  non-zero  value  of  meridional  flow  speed. 

Figure  1  illustrates  some  of  our  procedures  graphically.  The 
upper  left  frame  shows  a  Kltt  Peak  magnetogram  on  December  30,  1976. 

A  large  bipolar  magnetic  region  (^R)  Is  visible  in  the  southern  hemis¬ 
phere  on  this  otherwise  quiet  day  near  sunspot  minimum.  The  upper 
middle  frame  shows  these  fields  after  suppressing  the  short  wavelength 
components  of  their  Fourier  representation.  Obviously  on  this  day 
the  Kltt  Peak  zero  level  was  biased  in  the  negative  (blue)  direction 
and  the  positive-polarity  north  polar  field  is  not  visible.  The  upper 
right  frame ,  displaying  the  normal  component  rather  than  the  llne-of- 
slght  component  of  Held,  shows  the  same  BMR  deposited  within  the 
artificial  background  field.  Relatively  few  BMRs  have  contributed 
to  this  background  field  at  this  early  stage  of  the  new  sunspot  cycle. 

The  lower  frames  in  Figure  1  show  the  situation  28  days  later 
on  January  27,  1977.  The  lower  left  and  center  frames  show  the  raw 
and  smoothed  Kltt  Peak  measurements ,  respectively.  The  lower  right 
frame  shows  the  BMR  and  its  surroundings  as  they  would  have  evolved 
with  a  diffusion  constant  of  680  km^/sec.  At  this  time ,  a  newly 
emerged  BMR,  visible  in  the  magnetogram' s  northern  hemisphere,  has 


noc  yet  boon  deposited  in  the  simulator.  A  recently  deposited  doublet 
Is  visible  west  of  the  large  BMR  in  the  southern  hemisphere.  Our 
measurements  are  based  on  a  comparison  of  the  flux  distributions 
only  in  the  immediate  vicinity  of  the  evolving  BMR  and  not  in  the 
surrounding  areas  where  disagreement  is  both  obvious  and  expected. 


Fig.  1;  Observed  BMR  (upper  left  and  center)  is  deposited  Into  simu¬ 
lated  field  (upper  right)  ,  evolved  (lower  right)  ,  and  compared  with 
later  observed  field  (lower  left  and  center). 

In  conclusion,  we  should  like  to  emphasize  that  our  study  is  just 
reaching  the  point  where  most  of  the  technical  problems  have  been 
overcome.  While  our  results  to  date  are  tentative,  they  are  encour¬ 
aging  and  in  particular  give  a  diffusion  constant  that  is  consis¬ 
tent  with  Leighton’s  original  estimate. 

We  are  grateful  to  C.  Keppel  for  technical  support  ,  and  to  i.W. 
Harvey  for  Kitt  Peak  magnetograms.  This  work  was  supported  bv  NASA 
(DPRs  W  14,  429  and  W  14,777),  and  by  AFGL  (FY7121-82-03062). 
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